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Abstract: We compute the fully renormalized one-loop effective action for two interact- 
ing and self-interacting scalar fields in FRW space-time. We then derive and solve the 
quantum corrected equations of motion both for fields that dominate the energy density 
(such as an inflaton) and fields that do not (such as a subdominant curvaton). In partic- 
ular, we introduce quantum corrected Friedmann equations that determine the evolution 
of the scale factor. We find that in general, gravitational corrections are negligible for the 
field dynamics. For the curvaton- type fields this leaves only the effect of the flat-space 
Coleman- Weinberg-type effective potential, and we find that these can be significant. For 
the inflaton case, both the corrections to the potential and the Friedmann equations can 
lead to behaviour very different from the classical evolution. Even to the point that infla- 
tion, although present at tree level, can be absent at one-loop order. 
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1 Introduction 

High-energy processes in the early Universe must ultimately be described in terms of inter- 
acting quantum fields evolving in a curved space-time background. Often the dynamics is 
described by classical field equations of motion and the Friedmann equations for the scale 
factor, which are solved analytically or numerically. A plethora of models exist, where 
scalar field potentials are designed to generate a certain behaviour and certain physical 
phenomena. Although rarely stated explicitly, the understanding is that the classical field 
is really the homogeneous one-point function, or mean field, of a quantum field operator, 
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rolling in a quantum effective potential. This connection is often left unclear, and it is not 
obvious that for a given effective potential, there exists such an underlying renormalizable 
theory. In fact, consistent renormalization is rarely addressed in this context, opening up 
a number of pitfalls. 

In recent years, interacting quantum fields have become the subject of renewed interest 
in the general context of inflation, with the advent of searches for non-gaussianity in the 
Cosmic Microwave Background, curvaton models [1, 2], preheating [3] and dissipative in- 
flation models [4]. A number of approaches and approximations to the quantum dynamics 
have been employed, often motivated by practical tractability and features of the specific 
problem allowing for certain simplifications. In many cases, quantum corrections of one 
type or another are added to the classical dynamics, while others are neglected. It is not 
always clear that such choices are systematic, i.e. in terms of an improvable sequence of 
approximations . 

Quantum field theory in curved space-time has a long history (see [5, 6] and references 
therein) and since it is a first principle method, the difficulties in cosmological pertur- 
bation theory related to renormalization [7] and choosing the correct sized box for the 
curvature perturbations [8], are not present. Current work is focused on applications of 
the CTP ("closed-time-path") or Schwinger-Keldysh ("in-in") formalism to such spaces. 
In particular, the 2PI ("two-particle irreducible") effective action formalism was pioneered 
in this context by Calzetta and Hu [9] (see also [10]). This is a very powerful framework, 
where evolution equations for the mean field and the propagator are solved self-consistently, 
resumming a large class of perturbative diagrams. 

A drawback is that beyond leading order (1/N-LO or Hartree), the 2PI equations 
involve a memory integral, in addition to the doubling of time-labels of the propagator. 
This makes them numerically very hard to solve. Apart from [11], work has therefore 
concentrated on the leading order approximations, often in strict deSitter space [12-14] but 
also for self-consistent FRW Universes following from the Friedmann equation [10, 15, 16] 
(see also [17, 18]). These approximations are Gaussian in the sense that the only non-zero 
connected correlators are the mean field and the propagator 1 . Gaussian dynamics possesses 
a non-physical fixed point, and does not allow for dissipation and thermalization (see for 
instance [19]). All of these comments apply to Minkowski space applications as well, and 
are not specific to curved space-time. 

The issue of renormalization in this context has also received abundant attention. The 
benchmark method is adiabatic regularization [20], see also [5, 6], where counterterms 
are computed typically at the level of the evolution equations for the interacting fields; 
but also to renormalize the quantum expectation value of the energy-momentum tensor 
in the semi-classical Friedmann equation [10]. (For a different renormalization method 
functioning at the level of the evolution equations see [15, 16]). Adiabatic regularization is 
a very elegant and intuitive approach, whereby computing the adiabatic vacuum to fourth 
order in derivatives of the scale factor, all divergences can be identified and subtracted. 
One must then have some prescription to fix the finite parts of the counterterms. 

x It is however possible to obtain non-Gaussian curvature perturbations from Gaussian field perturbations. 
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Renormalization of 2PI equations beyond LO is notoriously difficult [21], and it is 
not obvious how to directly apply adiabatic regularization. This may be circumvented in 
deSitter space, where the Friedmann equation becomes irrelevant, and where we know the 
vacuum wherein to compute the full counterterms [14] . But in general, the lack of a simple 
renormalization procedure restricts the application of the 2PI formalism to cosmology. This 
needs to be resolved in future. 

However, it is possible to side-step some of these issues, if one is primarily interested 
in the quantum corrections to the mean field. This applies to the case where the field 
in question evolves in a given FRW background (we shall refer to this as a spectator or 
curvaton field), and to the case where the field itself dominates the energy density (we 
call this a participating or inflaton field), and where we therefore need to solve a quantum 
corrected Friedmann equation in addition to the mean field equation of motion. For this 
purpose, it is sufficient to consider the 1PI effective action. 

In this paper, we compute the fully renormalized one-loop 1PI effective action for a 
two-scalar theory 2 . This allows us to study inflation and curvaton evolution coupled to 
another field, and the effect of quantum fluctuations on the dynamics. As an alternative 
to adiabatic regularization, which was used for similar problems in [22, 23], we opt to 
renormalize at the level of the action using an expansion of the heat kernel [24-26]. This 
has the advantage that we can impose a set of renormalization conditions, which allow us 
to fix the counterterms completely. From the renormalized (and therefore finite) effective 
action, we derive evolution equations for both the mean fields and the metric by variation. 
These we solve numerically. 

The paper is organized as follows: In section 2 we write down the effective action. At 
this point we will simply state the result, which in the most general version includes inter- 
actions between the fields as well as self-interactions, and allows for non-minimal coupling 
to gravity. To illustrate the main effects of quantum corrections, we then immediately spe- 
cialize to two simpler cases for which we solve the evolution equations in section 3; first a 
curvaton in a quadratic potential coupled to another field in its vacuum; then an quadratic 
inflaton also coupled to such a field. We do not restrict ourselves to deSitter space. 

In section 4 we present in complete detail how we arrived at the effective action stated 
in section 2 including the crucial issue of renormalization, as well as a discussion of the 
approximations made and the limitation that apply to the result. Further technical points 
can be found in the appendices. We conclude in section 5. 

2 Effective action of two-scalar field model 

We consider two real scalar fields a and 4> in what will eventually be a flat FRW background, 
but which for the moment is described by a general metric g^ v '. Depending on the context, 
the fields could be the inflaton, the curvaton or some other scalar fields, but at this point 

2 This approximation is formally equivalent to the "less than leading order" truncation of the 2PI effective 
action, i.e. where no 2PI diagrams are included in the functional IV At this level of truncation, the Keldysh- 
Schwinger contour is less crucial, but at higher loop order this will come into play. In the language of [9], 
the propagator G is G ++ . 
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they are general and on equal footing. Our convention for the metric is (—,+,+,+) and 
the Riemann and Einstein tensors 3 are defined according to the " + " conventions of [27]. 
The bare (classical) action reads 



S[<i>,a,g^]=S m [<i>,a,g^] + S g [g^) 



(2.1) 



where 



S m [<P,a, 9 n = j d 4 x^ 
1 



1 Tfl± 1 /? 

-g^d^dA + v^ 2 - -^4? - -(1 - £ )-</> 2 



g(f) 2 cr 2 X a cr A 



4! 



4! 



(2.2) 



and 



S 9 [ 9 H = J d A x^-g 



A + aR + /3R 2 + eiC 2 + e 2 G + kDR 



(2.3) 



We have introduced the square of the Weyl tensor C 2 and G is the Gauss-Bonnet density, 
defined as 



C 2 



- R 2 — 2 R pv f? _i_ w vp(7 R 



G = R — AR^R^ + R^ pa R 



(2.4) 



in terms of the Ricci scalar R, the Ricci tensor R pu and the Riemann tensor R pV(XT . We also 
have the operator □ = V^V 1 = {gl^ 1 ^ 2 d^dgl 1 ^ 2 d p ) . We assume that m 2 , m 2 , A^, A CT and 
g are all positive. Negative masses squared could be considered, leading to spontaneous 
symmetry breaking, which is especially interesting in curved space [28], but we will refrain 
from doing so here. The couplings r/ and k multiply total divergences and in unbounded 
spaces the boundary terms vanish. We also used the parametrization for the £ for which 
£ = or 1 means the conformal or minimal coupling respectively. /3 = ei=e2 = K = and 
a = I/WttGn (Gn is the Newton constant) is the Einstein-Hilbert action with two scalar 
fields. In terms of the "usual" cosmological constant A', we have A' = SttGnA. To each 
of the 15 parameters corresponds a counterterm, which will come into play below, when 
we discuss renormalization. 

The 1PI effective action F[<p,a, g^] is a functional of the mean fields 4 



a(t) = (a(x,t)), <P(t) = {<P(x,t)), 



(2.5) 



which for a homogeneous state have only time dependence. We will treat g^ u as a clas- 
sical external field. There are therefore no quantum gravity corrections included in our 
treatment, only corrections from loops of the scalar fields, which as we will see include 



pM I pM rff 

1 a/3, 7 ~T L ct/3 1 - 7° 



4 We will use the same symbols a 



2 ^9^^ ^/Af- 

' to denote the mean fields. It will be clear from the context whether 



pm r" n 

^ (77 /3ai 



we speak about the mean fields or the field operators a(x,t), <j>(x,t). 
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gravitational operators. Once the effective action is known, it is straightforward to derive 
the equations of motion by variation 



Sm<r,g^] 



8g 



fJ,V 



ST[^a,gn _ Q ST[J>,a,gn 



8a 



0. 



(2.6) 



The first of these amounts to a quantum corrected Einstein equation, and the last two are 
quantum corrected field equations of motion. 

We compute V at one loop order, and fix the counterterms in such a way that at the 
point g^ u = rf^ , a = (f> = the renormalized couplings are equal to the tree- level ones. 
We postpone the details of the computation and renormalization of the effective action 
to section 4. Truncating a derivative expansion in the way described there, we have the 
general result 



>(0) , r W 



T[ct>,a,g^} = J d 4 x^g-£ eff = J d A x^ 
with a part that looks like the tree-level action (but now in terms of mean fields) 

+A + aR + (3R 2 + Ge 2 , 



AO) 

^eff 











if 





4! 



A„a 4 



b 2 a 2 



4! 



C 



and a quantum correction 

S/ = 64^ { ki I 69 {9 + 3 (A<T + X<t>)) ^ + ^ + R2 ~ 4 (^ m - + H 
+9 (g 2 + Aj) 4 + 12 (gml + m 2 A^) 4> 2 - 6 (<?£, + A^) Rcf 2 

+9 (g 2 + A 2 ) a 4 + 12 (grn 2 + m 2 A CT ) a 2 - 6 (^ + A CT £ CT ) ^ 



+ 



+ 



G 



g 2 a 2 4> 2 ml 



1 



360 4(m 2 - ml) 2 
G 



- + 



+ 



+ 



g 2 a 2 4> 2 m 2 a 



+llo 
12 



360 4(m 2 — 

M 2 



1 / 2 , ° X « , 

- m„ H 

2 V 2 2 





log 


(M 2 _M\ 








log 


(MlM 2 




V "4 



3a 2 ( g + A CT ) + 3</> 2 (<? + A*) + 67T7 2 + 6m^ - Rfa + ^)] (M 



(2.7) 
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We have defined 



M. 
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where the mass-eigenstate fields are given by 



1 



<s>l = ^g(* 2 + <P 2 ) + - , 



+ X a a 2 



+ y nr, - hi 
1 



1 \ : 



+- a 4 (g - KY + ^(g- Xcf,) 2 + 2^V {g 2 + (A CT + X^) g - X a X 



1/2 



(2.11) 



One notices that the effective action (2.8), (2.9) has a very simple structure. As for the 
classical action, it is completely symmetric in <p and c, and in particular the second and 
third lines of (2.9) are the same with o a, as are the next two. 

We have specialised to FRW space-time, defined in terms of the scale factor a(t), 



ds z 



-dt 2 + a 2 (t)dx 2 , 



(2.12) 



where we can define the Hubble rate H = a/a and the higher order gravitational tensors 
reduce to 



C 2 = 0, 



G = 24 



i 2 •• 

a \ a 



R = 6 



a / a 



+ 



(2.13) 



As explained in section 4.3, partial integration of □ x field type terms multiplying the 
logarithms in (2.9) reveals that they are beyond our level of approximation 5 . Also, the 
□ x field terms that would appear in (2.8) integrate to zero. These have all been discarded 
already along with C 2 . For the full general results, see Eq. (4.31) and (4.32). 



Gravitational corrections enter through the second order operators R 



and the 



fourth order operators G, C 2 , R 2 , HR 



7^- 



We have truncated keeping this order, which 



we refer to as 0(R 2 ) (see section 4). Ignoring these, the expression reduces to a Coleman- 
Weinberg-type effective action. In particular when also setting g = 0, it is easy to see that 
we get the sum of two decoupled Coleman- Weinberg potentials for a and <fi. The quantum 
correction C^) f is suppressed by the phase space factor l/(647r 2 ) and is in this sense small. 



-eff 

Also reinstating physical units, h would appear in front of C 



(i) 
eff 



2.1 A specific case: FRW, minimal coupling to gravity and no self-interaction 

The above expression in all generality may be applied to a number of cosmological scenarios 
and phenomena, but in the following we will simplify things somewhat by taking the values 



£<t = = i) P — e i — £ 2 — x a — x^ — o, 



(2.14) 



3 This is true in unbounded space. Otherwise, there may be additional boundary terms. 
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and use the symmetry of the potential to solve for one of the fields in its equilibrium 



= 0, (2.15) 

leaving the other field a to evolve freely. We emphasize that this is simply an appropriate 
parameter choice for the particular application we have in mind here, and not a further 
approximation. This allows us to write 

£eff = — — 2 ^a 2 + A + aR 



+ 



2 R + 9<? 2 \\ G 
m ^"6 + ~J + 180 



log 



m <j> 6 + 2 
o 

mi 



2 



2 R\ G 

mi H 

6 J 180 



m 2 - ^ 



log Mbr^ ■ (2-16) 



In fact, the entire last line can be discarded since an expansion in R/M± reveals that it is 
also 0(R 3 ), with the understanding that the logarithm must be finite, i.e. i?/(6m 2 ) < 1. 
Similarly from the first logarithm, and anticipating that a is to oscillate through zero, we 
must also demand (R/6mV) < 1. We will see in section 4 that this is partly a result of our 
choice of renormalization scale. We also mention here that for the derivative expansion to 
be reliable, one must require <? < 1. 

3 Quantum corrected field dynamics 
3.1 Dynamics of spectator fields 

We start by considering the case where the fields a and (ft do not dominate the energy 
density. The obvious example of this is the subdominant curvaton during inflation. We 
assume that there is some other energy component determining the evolution of the back- 
ground metric, in FRW the scale factor a(t). We will consider three typical cases, radiation 
dominated, matter dominated and de Sitter space, where the scale factor evolves as 



OradW = o ( — > cWt(i) = a ( — — -J , a dS {t) = a e Hot , (3.1) 



and we have the Hubble rate 



and the higher order operators 



-Rrad — 0, -Rmat — _|_ ^2 ^ dS ~~ ^Hq, (3.3) 
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Figure 1: The evolution of the spectator field a(t) for different values of do and different 
approximations, in a matter dominated Universe. We use m^/m^ = 2, g = 1, and Ho/m a = 
a/1/2. The inset shows the same thing but rescaled by the initial values ctq- 



and 

3 64 
G ™ d = ~W+W' Gmat = "27(t + t ) 4 ' GdS = 24 ^o- ( 3 - 4 ) 

Since the metric is given and we have initialized (f> in its minimum, eft = 0, we only have to 
consider the equation of motion for a which reads 

1 ( i G (m 2 R i 9" 2 \ 2 \ 

•• , qtt- , 2 1 I / 2 n I n 2 \ i V180 V'V 6 2 ) ) 

a + Wa + m a o = -^l-{2m 4> -R + Zgo)+go 2 _ R g 

k "V 6+2 

" ¥ K " e + — ) lQ g ( ^2 J j = 0' (3-5) 

At tree level this reduces to 

a + 3H& + m 2 a a = 0. (3.6) 

and there is then no dependence on g or m^, and a does not feel the presence of the 0, to 
which it is coupled. 

We solve the evolution equations in the three cosmological backgrounds, identifying 
four levels of approximation: Tree level is (3.6), order H° ignores gravitational operators 
in the quantum corrections, order H 2 includes all occurrences of R and finally order H 
also includes G. We choose the parameters m^/m^ = 2, Ho/m^ = y/l/2, g = 1- m a sets 
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Figure 2: The evolution of the spectator field o{t), for different values of do and different 
approximations, in a radiation dominated Universe. We use m^/m^ = 2, g = 1, and 
Ho/rrifj = a/1/2. The inset shows the same thing but rescaled by the initial values <tq. 



the scale, and Hq is chosen so that the argument of the logarithm is always positive. In 
particular Ro/(6mi) = 1/2. 

We then vary the initial amplitude ao/m a , and since (3.6) is linear the evolution is 
simply rescaled. The full dynamics is no longer linear, and such a rescaling no longer works. 

In Fig. 1, we show the evolution for a matter dominated background, in Fig. 2 for 
radiation and in Fig. 3 in de Sitter space. We see that the field performs a damped 
oscillation, which in the de Sitter case is even over-damped. Choosing ao/m a = 1, 10, 20, 30, 
we get the black red, green and blue curves, respectively and where the dashed line signifies 
the classical tree level result and the full lines the quantum corrected ones. In the insets, 
we show the same curves rescaled with <to, to demonstrate the deviation from the linear 
evolution. 

The quantum corrections are numerically relatively small at these values of parameters, 
which is simply due to the overall factor l/(64-7r 2 ). We also find that different approxima- 
tions of order H 0,2 ' 4 as defined above are essentially identical (they all fit within the full 
lines). This means that the quantum effects are due to the non-gravitational corrections, 
i.e. the Minkowski space effective potential. 

Our choice of renormalization conditions (see section 4.4) tunes the counterterms so 
that the effective potential matches the tree level one close to a = 0. But for field values far 
from zero, the effective potential is much steeper and non-quadratic, leading to non-linear 
oscillation. These non-linearities are a quantum effect. 

Choosing a different renormalization point changes the effective potential to match at 
some other value of a, but then that effective potential has large corrections near a = 0. 
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Figure 3: The evolution of the curvaton field cr(t), for different values of 00 and different 
approximations, in a de Sitter Universe. We use m^/m^ = 2, g = 1, and Ho/ma = y/l/2. 
The inset shows the same thing but rescaled by the initial values ctq. 



The different choices are physically distinct, and should be based on what the value of the 
physical, renormalized parameters are taken to be. 

This concludes our study of spectator fields in a given gravitational background. We 
now turn to the case where a and <p dominate the energy density. 

3.2 Quantum corrected Friedmann equations 

The quantum effects in the effective action also lead to corrections to the Friedmann equa- 
tions, which follow from taking the variation 6 

— l — c J = 0. 3.7) 
Sgf a> 

Since the effective action has curvature terms up to 0(R 2 ) the Friedmann equation is now 
fourth order in time derivatives (see 4.5). One may reduce the fourth order Friedmann 
equation to a second order one keeping all the physical solutions [29] , but the procedure is 
rather involved. Here we settle for 0{R). We will also restrict ourselves to the simplified 
case (2.16) for the field dynamics, but an analogous procedure of course applies to the 



'Note that one must vary w.r.t. the general metric g^ v and only afterwards specialize to FRW. 
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general result (2.8) and (2.9). We find from the "00" component 
■ 2 ( Gn 



12vr 



a( G N a 2 (,.9<T 

+ -\ TTj 9° log 1 + t— j 

a [ Yin a \ 2m A , 



A' 8wG N f& 2 ml 2 
3 3 1 2 2 a 



Gn 
12n 



ga 2 (Am 2 +3ga 2 ) i oa 2 2 / 5 <7 2 



16 



2m^ 



(3.8) 



This result has a straightforward interpretation. When neglecting the loop contribution 
(terms proportional to Gjv/(127r)), the equation reduces to the well-known classical Fried- 
mann equation 



d 2 8ttG n (a 2 ml 2 \ A' 



8itG n A' 



(3.9) 



where p = Too is the classical energy density in the fields. The quantum corrections involve 
quite non-trivial field-dependence and in addition have explicit dependence on H = a/a. 
We note that occurrences of R 2 (for instance coming from expanding the logarithms) have 
been neglected through neglecting derivative orders beyond two. 

Similarly, for the "ii " component we find, after using (3.8) to simplify the expression 



a] 24vr 



{2ml + ga 2 ) log 1 + 



9° 
2m? 



90 



h! 4ttG n 
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2d 2 + m 2 a 2 ) 



+ 



Gn 
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ga 2 Um\ + 3ga 2 ) 



9 2 a 2 cr 2 



+ 



a( G N a 2 / ' go 2 



2 , 9<? 



2\ 2 



[go 2 + 5 <7<j) ) log ( 1 + gn 



2m 2 
(3.10) 



Again, the interpretation is clear, and neglecting the quantum corrections, we find the 
standard expression 



d AttGat , 9 9 9 N A' 
- = — (-2. 2 + ^ 2 ) + y 



AttGn , „ A' 
■^-(p + 3 P ) + y , 



(3.11) 



with the classical pressure p = Taj a 2 (no summation). The expression now also involves 
second derivatives of the field. There is again a term proportional to a/ a, which vanishes 
in the classical limit. 

If we choose to separate the effective lagrangian (2.16) in the form R + aC 1 (a, R), 
Eqs. (3.8) and (3.10) correspond to the classical Friedmann equations, where we replace 
the right-hand side by the appropriate components of the quantum corrected energy- 
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Figure 4: The evolution of the inflation field a(t) rescaled by its initial value <ro, for 
different values of <7o and different approximations, when the scale factor is given by the tree 
level (full lines) and quantum corrected (dashed) Friedmann equation. We use m^/m^ = 2, 
<? = !■ 



momentum tensor 



2 5 f d 4 x^gC rest (a,R) 



T -- ~ — V- <3 ' 12) 

However, it is not obvious that this is the same as replacing the right-hand side of the clas- 
sical Friedmann equation by a renormalized quantum expectation value of the components 
of the energy momentum tensor (see for instance [11, 12]). For example, if one starts from 
the renormalized effective action, the covariant conservation of energy-momentum follows 
directly from the equation of motion for g^ u (2.6), if energy is defined as in (3.12). In 
contrast, if one starts from the classical Friedmann equation with the energy- momentum 
repalced with the quantum expectation value and then performs renormalization, the co- 
variant conservation becomes a non-trivial issue [22]. 

3.3 Self-consistent dynamics of field and scale factor 

In the following, we will solve the equation of motion for a together with (3.10) for the 
tree level case (3.6) and the quantum corrected equation (3.5). We take the parameters 
Mpi/m a = 100, g = 1, m^/m CT = 2 and vary ao/m a . We define the reduced Planck mass 
as Mpj = 1/(8it 2 Gn). The initial value of the Hubble rate Hq is given by (3.8), and we 
simply set do = 0. 

Fig. 4 shows the field a using different initial amplitudes ctq, and rescaled by this 
amplitude. Whereas for the spectator field above this rescaling is exact for the tree level 
equation, this is no longer the case here, since the Friedmann equation has a non-linear 
dependence on the initial amplitude. In the left-hand plot we show ao/m a = 1 and 10, and 
only for the larger amplitude is the quantum corrected evolution (dashed line) discernably 
different from the tree level evolution (full line). In the right-hand plot we see that this 
trend persists for even larger initial field amplitude. 

In Fig. 5 we show the corresponding scale factor a(t), for all the different initial am- 
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Figure 5: The evolution of the scale factor a(t) from the tree level (full lines) and quantum 
corrected (dashed) approximations. 




plitudes and for the tree level (full) and quantum corrected (dashed) equations. All curves 
turn out to evolve as matter domination oc i 2 / 3 , but with different rates. Again the quan- 
tum correction is negligible for the smallest amplitude (black lines), whereas there is a 
significant effect at larger amplitudes. 

3.4 Inflation 

Of particular importance for cosmology is the possibility of inflation at the very earliest 
stages of the evolution of the Universe. In a simple slow-roll computation, we have for the 
tree level quadratic potential 

1 2 2 2M P1 , s 

V = -mla 2 , 6 = — (3.13) 

Hence, neglecting quantum corrections, inflation will arise for do > v2M p i, when e < 1. In 
Fig. 6 we show the evolution of a(t) for o"o = 10M p i for early times, for M p \/m a = 100 (left) 
and 10 5 (right). Indeed in the tree-level case (or non-interacting, g = 0, red line), inflation 
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is prolonged and leads to exponential growth of the Universe. Now we gradually add the 
interactions by increasing g (full, black lines). We see that inflation becomes weaker and 
weaker, and finally disappears for g ~ 1 (left) and g ~ (right). Hence because of 

quantum effects, interactions with another field (ft may prevent inflation from taking place, 
at least in this case for oq < 10Af p i. This is a well-known result, which has now also been 
confirmed in the present framework. 

The value M p i/m a ~ 10 5 is what is required for the spectrum of the CMB to be 
correctly reproduced by quadratic inflation (see for instance [30] ) . If the inflaton is coupled 
to other fields, this conclusion must be re-evaluated in the context of quantum corrections 
to the potential. Since an analytic solution of the slow-roll field evolution from the end 
of inflation back to the era of horizon crossing is not readily available, we will postpone a 
detailed analysis of this point to future work. 

3.5 Field fluctuations from the propagator 

A crucial observable for inflationary physics is the spectrum of density perturbations, which 
arise from the quantum fluctuations of the inflaton or a spectator field like the curvaton. 
From the effective action point of view these are encoded in the field propagator (or two- 
point function), which follows from taking the second variation of T w.r.t. the fields, i.e. 

^ = 6xi5x2 . (3-14) 
where xi i can be (ft or a. This gives a matrix equation for the propagator 

^ 1 ^ Gaa Gfjcf, \ , / //x 5(x — x") 



g-Hx,x')g(x', x ") = %° (x, x >) I ;r I = ; > (3.15) 



where an integral over the x' indexes is left implicit. At the one-loop level employed here, 
the propagator equation is equivalent to solving the equation of motion for each field mode 
fk(t) of the mass-eigenstate fields 

/d 3 k r i 
[4fk(t) + (4) t (/ fe ± )*(*)j . (3-16) 

where at are creation operators and the equation of motion is 



„2 

r2. 



fk ~ 3Hf k - \h - M 2 a {a, <t>)f k = 0. (3.17) 



for some M%(a,(p). For sufficiently large momenta k the mass term can be neglected as 
is usually done, but at horizon crossing H = k we have H > M(cr, (ft), and so formally 
the truncation leading to the effective action in (2.8, 2.9) no longer holds. But since 
the equation for the propagator in that limit is the free field equation used in standard 
computations of the spectrum we just conclude that we have no quantum corrections for 
this observable, other than what enters via the altered evolution of a(t) and cr(t). 
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At higher loop order, additional interactions will make the propagator equation non- 
trivial. One way to introduce these is through the 2PI formalism. For the 1PI effective 
action, one first solves for the mean field, and this is then put into the propagator equation 
which is then solved. For the 2PI formalism, both mean field and propagator (and scale 
factor) are solve for simultaneously, and the equations are coupled. This provides the self- 
consistent solution which is formally a resummation of a particular sub-set of perturbative 
diagrams. As we discussed in the introduction, technical hurdles mean that the 2PI for- 
malism has not been completely adapted to expanding backgrounds beyond LO (although 
see [11]). 



4 Computation of the effective action and equations of motion 
4.1 Effective action for <p 4 theory 

We will first discuss the relevant points in deriving the curved space effective action for (p 4 
theory to one-loop order. Our discussion follows closely that presented in [25]. For more 
details, see also [5, 24, 26]. We start from an action with a matter part 



Smi^gH = J d 4 x^~g 
and a gravitational part 



1 



R 



A 



-g^d^d u ip + ryfV - - ~(1 - £)^ 2 - ^ 



(4.1) 



S g [ 9 n = / d 4 xV=g~ 



A + aR + PR 2 + eiC 2 + e 2 G + kDR 



(4.2) 



Our aim is to write the effective action as a loop expansion, 

T[<p,gr] = J d 4 x^C eff =T^[ l p,g^]+T^[ V ,g^} + --- , (4.3) 

where the metric g^ v is considered as a classical field. The standard method is to start 
from the generating functional for connected diagrams, use an expansion of the action to 
quadratic order, discard the one-particle tadpoles and finally perform a Legendre transfor- 
mation to obtain the results [31] 

r^foflH = S m [<p,<r] + S g [ 9 n + SS[ V ,gn, rW^sH = ^Tr log G~\x,x'), (4.4) 



where G{x,x') is the Feynman propagator satisfying 



1 



21 



n + m z + -R(i-0 + ^Y 



G(x, x) 



5(x — x') 



(4.5) 



One may choose different boundary conditions (i.e choose different vacuum states) for the 
propagator in (4.4) but this will turn out irrelevant for our approximation. The counter 
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terms enter as, 

5s[<p,<r]-- 



d x-y/— ~g 



$Z o <- r-, 2 2 c^-R 2 ^ 4 

—g^d^dvtp + <5ryD(/7 - — c/? 2 + <5£— y> 2 - — 



5A + 5ai? + <5/3i? 2 + <5eiC 2 + 5e 2 G + 5k\JR 



(4.6) 



and are introduced for the purposes of renormalization. In order to find an expression for 
we use the heat kernel method introduced for curved spacetime in [32] and write the 
trace of a logarithm in (4.4) as a proper-time integral over a yet undefined kernel function 
K 

■ l ~ n I d n x^ J ^-K(t;x,x). (4.7) 

Because of the divergent behaviour that occurs in four dimensions for , we have dimen- 
sionally regularized the above integral to the dimension n = 4 — e and added an arbitrary 
scale [i in order to maintain the proper dimension of the action. In appendix A we present 
the heat kernel solution in full detail. Here we simply state the result which reads 



-Trlog G~ 1 (x,x') 



K(t; x, x) = i 



f2(r; x, x)e 



-iM 2 T 



where M is the effective mass 



M 



m 



(4.9) 



and the O's have a small proper-time expansion 



0(t; x, x) = ak(x, x)(ri) 
fc=o 



(4.10) 



Substantial literature exists on these coefficients [33] and [34] and results for the first four 
can be found for example in [35]. Inserting the expansion (4.10) into (4.8) and performing 
the proper time integral in (4.7) we get for the effective action 



d n x, 



1 



n-4 oo 



2(4vr) n /2 ^ n J 



Y,M 4 - 2k a k (x,x)T(k-n/2). (4.11) 



k=0 



Since we made use of a small proper time expansion to approximate the entire integral 
over the heat kernel, some comments on the validity of the procedure are in order. As is 
evident from the divergences of the Gamma function on the right hand side of (4.11), the 
small r part of the heat kernel probes only the local region of space-time, and may miss 
some global aspects. For example, as is argued in Appendix A boundary conditions are not 
seen by the expansion (4.10). The result for the full heat kernel is in fact non-local and has 
been developed recently, [36] and references therein, and applied to quantum scalar fields in 
[37]. The ansatz (4.8) is, however, a valid approximation for the full effective action when 
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the fields are slowly varying [38], [39] and [40] 7 . The coefficients 8 a& consist of geometric 
tensors and derivatives of ip 2 and from (A. 13) can be read to be ao = 1, ai = and 



/ \ a ,— , 9 - \ n R 3C 2 — G 

a 2 (x,x) = Dip + (5f + 1) 1 , 

12 ^ v s '180 360 



(4.12) 



where we have defined C 2 and G in (2.4) 9 . The approximation to be used in this paper 
is to include only the first three terms of the right hand side of (4.11) in rM, i.e. ao, a\ 
and d2, neglecting a% and beyond, and then subsequently neglecting terms in a 2 that also 
appear in 03 (see section 4.3). This gives us the effective action for a (/? 4 theory to one-loop 
order 



r>, <n = s m [<p,<r] + s g [ 9 n + 5s^ g n 



+ / d U Xy 



64vr 2 







3" 


| — M A 


Hf) 








~ 2 



,'M 2 
2a 2 (x,x)log ( 



where we defined a new arbitrary scale \x 



log(/i 2 



n 



7 + log(47T// 2 ). 



(4.13) 



(4.14) 



From the right hand side of equation (4.13) one can immediately see that the (regularized, 
but still unrenormalized) effective action consists of the well-known Coleman- Weinberg 
effective potential [41] and corrections coming from field derivatives and gravity contribu- 
tions. We reserve the closer investigation of the implications of our truncation in (4.13) to 
section 4.3. Renormalization will be discussed in section 4.4. 

4.2 Effective action for the two-scalar field model 

Guided by the results for the ip A model, we can straightforwardly derive the one-loop 
effective action for our two-field theory in (2.2). The one-loop contribution will be the sum 
of two functional determinants 



T^[cj>,a,g^} = *-Tr log A+(x,x') + ^Tr log A_(x,x'), 



(4.15) 



where A± are the eigenvalues of the mass matrix 



M 



A^ 2 



□ + m} + (l-fc)£ + '=- + ^ 
g 



kcrc 



□ + m 2 + (l_£ T )# + flf + 



A CT a 2 



(4.16) 



7 Fields are smooth on the scale of M 2 . 

8 Not to be confused with creation operators au for quantum fields. The have adopted established 
notational conventions. 

9 Even when we have analytically continued to n dimensions, we continue to use the four dimensional 
definition for C 2 . 
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We parametrize the eigenvalues as 



A+ = 



□ + 



mi + m 2 



2 2 



5(x — x') 



(4.17) 



where &± is defined in (2.11). Retracing the steps that led to (4.13) allows us to write the 
full effective action for two scalar fields to one-loop order as 



Y[<t>,a,g^] = S m [4>,a,g f11 '] + S g [g^ u ] + 5S[<p,cr,g^ u ] 



+ J d n xy 
+ f d n xy 



1 



64vr 2 
1 

64vr 2 



M% 
Ml 



log 
log 



M% 
Ml 



2a^(x, x) log 
2a^ (x, x) log 



Ml 
Ml 



with the natural definitions 



Mi 



m 2 , + ml 



2 ± 



af(x, x) 



1 

12 



+ 5 



12 



+ 1 



OR 3C 2 - G 



180 



+ 



360 



(4.18) 

(4.19) 
(4.20) 



From the expression for the eigenvalues &± in (2.11) it can be determined that at the limit 
g — > the result in (4.18) reduces to two decoupled effective actions of the form (4.13). 
At this stage our expression is still divergent in the limit n — > 4 since we have not yet 
renormalized, only regularized (with the scale jl still present). However as promised after 
(4.14), we first make some comments on the level of our approximation. 

4.3 Truncation of the effective action 

In our derivation for the effective action (4.18) we truncated under the assumption 

a 3 (x,x) 



Ml 



<C a,2(x,x) 



(4.21) 



where the terms in can be read e.g. from [35]. Schematically, for a c/? 4 theory, the 
neglected terms are of the type 



R 



U 2 R 



M 2 ' M 2 
It is also worth noting that terms of type 

' M 2 



M 2 



V^VHog 



m' 



RD log 



M 2 



M 2 
m 2 



DR 
M 2 ' 

R*_ 

M 2 ' 



(4.22) 



(4.23) 
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appear in the calculation (for instance from 02 after partial integration) and are of the 
same order as those in (4.22). These are also neglected. 

For our two-field model (4.18), the following terms are similarly neglected 



Ml ' 



± Ml ' 



R 



UR 

Ml' 



R 



DR 



"Ml"' 



□ 2 i? 



Ml' 



(4.24) 



For the specific choices ^ = £ CT = 1, = A<y = 0, we have 



10 



ml - m 2 , + . 



2 2,2 



(4.25) 



and 



M 



± 



6 2 



2 R , 
^ " 6 + ~T 



(4.26) 



From these considerations we can deduce that our approximations are valid as long our 
fields (a, <p and g^ v ) are slowly varying in agreement with our motivation for the anzats 
(4.8). 



4.4 Renormalization 

The renormalization procedure for the effective action is well-established, for example in 
the original paper [41], and for curved space [24, 25]. First, we briefly present the renor- 
malization procedure for the (p 4 theory in unbounded Minkowski space-time for a constant 
field. The effective action in (4.13) becomes 



J d 4 x^J^g C e f f = - J d 4 x^f^g V(<p) 



d 4 X\J —g 



m 



—p 2 - -p 4 + A 

2 Y AV 2 



5m 2 2 5X 4 , 
-f ~ -77 f + SA 



4! 



64vr 2 



m + X 



2\ 2 



log 



m 2 + Af 



(4.27) 



The counter terms 5m 2 , 5X and 5 A are fixed by a choice of renormalization conditions, for 
which we take 



m 



d 2 V{ip) 



d^ 2 



x 



d 4 V(p) 



<p=<pi 



dp 4 



-A = V(p) 



<P=<f3 



(4.28) 



where the fields <p n are arbitrary constant renormalization points. If we choose p\ = P2 - 
(fs = and further that A = 0, we get the standard result for the effective potential [42] 



-.r/ \ m 2 2 A 4 1 



2 *P 
m +A— ) log 



m 



+ Af\ Xip 2 (4m 2 + 3X<p 2 



m' 



(4.29) 



10 Actually, the eigenvalues $± in (2.11) acquire an absolute value coming from the square root of a 
square, but its sign will never become explicit in the effective action (2.9). 
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In a similar fashion, we can renormalize our two field effective action (4.18) in a general 
spacetime with the following renormalization conditions 11 



-ml = 


d 2 £eff 
Off) 


-ml 


_ d2 ^ff 




acr* 


—^<f> = 


d^eff 


-9 


_ d 4 C eff 


i-U _ 






deft 2 a 2 


6 


d(ft 2 8R 






A 


= C eff, 


a = 


d £eff 


6 


da 2 dW 


OR ' 


2/3 = 


8R 2 ' 


ei 


dC -ef{ 
8C 2 ' 


£2 = 


dG ' 


V<f> = 


dUcft 2 ' 




Meff 
dUa 2 ' 


fC = 


0or ' 



(4.30) 

where it is understood that for every renormalized constant the renormalization point for 
g^iu - g an j £ Qr ^ anc j ^ - g c ] losen t Q be zero. Furthermore, we will not have to deal 
with the infrared singularity mentioned in [41] resulting from the zero scale renormalization 
of the four-point functions, since we are not interested in the massless limit. Also besides 
renormalization, the massless limit would pose additional complications in curved space 
due to the lack of an exponential damping factor in the heat kernel anzats (4.8). The 
massless limit would probably require the use of the non-local heat kernel expansion as 
pointed out in [38]. One should therefore be very careful when taking the small mass 
limit of (4.29) (as was done in [43]), or when imposing renormalization conditions. Results 
for the counter terms defined in (4.30) can be found in Appendix B and by inserting the 
expressions into (4.18) and taking the limit n — > 4 we get the full results for the effective 
action 

4?/ = -\r 0,4^4 + - '4<p 2 - |(i - wf -p 2 - 

+A + aR + PR 2 + C 2 e 1 + Ge 2 + kDR, (4.31) 



11 There is no need for a condition for 5Z since the one-loop contribution has no momentum dependence. 
This will appear at higher loop order. 
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and 
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6g (g + 3 (X a + A*)) <rV + (£ + # 2 - 4 (£ CT m 2 + 



W = 1 f 1 

e ^ 64vr 2 \ 24 
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&7 2 ( ff + A CT ) + 3<£ 2 ( 5 + A ) 
2 fl(^ + ^)-2Dl(M+-M_; 



(4.32) 



If we are in an unbounded space, we can use partial integration to remove the terms with 
a □ since they are either total divergences or beyond our approximation, cf. section 4.3. 

4.5 The quantum Priedmann equation to order R 2 

Here we briefly go through the steps necessary for deriving the quantum Friedmann equa- 
tion used in section 3.2. The calculation of the Einstein equation 



5T[g 



(IV 



Sg 



(IV 



Sg 



(w 



-gc 



eff 



0. 



(4.33) 



for the action (2.16) is straightforward. Because of our truncation of the series in (4.11), 
there is no loss of accuracy if we first expand to 0(R 2 ) and then calculate the variation 
w.r.t. g^ v by using the formulae in appendix C discarding higher-order terms according to 
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section 4.3. This gives 
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(4.34) 



Using the tensor expressions from appendix D we get the Friedmann equations for the "00" 
component 
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(4.35) 



and for the "w" component 
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If we neglect the 0(R 2 ) contributions (four dots) and perform the redefinitions a — > 
l/(167rGjv) and A -4 — l/(87rG/v)A' equation (4.35) becomes (3.8) and if we further solve 
for a 2 /a 2 from (4.35) and substitute it to (4.36) we get (3.10). 

4.6 Equations of motion to order (R,g 2 ) 

As a final exercise, we also show the result of further expanding the equations of motion 
and the Friedmann equations (3.5), (3.8) and (3.10) to leading order in the coupling (still 
concentrating on our simplified model). The equations for a and a can then be written in 
a very physically intuitive form 

for a, and for a 



^fi+^4^+-f-^4^ 

a z \ 967rm^ / a\a2vKm^ ) 
a\ 96^771^ J a\o-48irmi ) 




where p and p are the tree-level energy-density and pressure defined in section (3.2). 

The leading order correction is proportional to cr 4 in all three equations (cr 3 for the 
cr-equation corresponding to a cr 4 correction to the potential) , but perhaps more surprising 
is that for the a equation it is also proportional to R, and hence vanishes in a radiation 
dominated background. On the other hand, we have seen in the previous sections that the 
quantum correction is insensitive to the gravitational operators. Hence the leading order 
term in g does not dominate, simply because R is in general small. We do not advocate 
the use of this approximation, but simply show it for illustration of the structure of the 
equations. 

5 Conclusion 

In the present paper, we have computed the renormalized one-loop 1PI effective action for 
two coupled and self-interacting scalars in a general FRW background. We used a gradient 
expansion for the effective action and truncated at 0(R 3 /M±), where M± are eigenvalues of 
the mass matrix. We renormalized at the level of the action, using a set of renormalization 
conditions to fix also the finite parts of the counterterms [24-26]. 

Subsequently specializing to a model of one field (cr) rolling under the influence of the 
quantum fluctuations of the other (0), we derived the quantum corrected field equations 
of motion. We also derived the quantum corrected Friedmann equations (3.8) and (3.10), 
which are in general different from what is commonly referred to as the "semi-classical" 
approach. 
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We solved the equations of motion for the case where a does not dominate the energy 
density and the evolution of the scale factor is given. We found that the quantum effects can 
be significant, but that gravitational operators can be largely neglected. This is important 
for instance for scenarios involving subdominant curvatons [44]. Our results improve on 
the discussion in [43] , and we have stressed how important it is to renormalize correctly to 
get the right physical result, using appropriate renormalization conditions. 

We then considered the more involved case, where an oscillating a dominates the energy 
density, and solved the coupled system of field equation and Friedmann equations to find 
the evolution of both a and a(t). We again found that the effect of including quantum 
corrections could be significant. This will be important for the reheating and preheating 
process after inflation. 

In a quadratic potential, a large initial value of a classically leads to slow-roll inflation. 
We found that including quantum corrections shortens the period of inflation, and may even 
prevent it from happening. This is just one instance of the general result that quantum 
corrections spoil the flatness of the inflaton potential. Such effects are less pronounced for 
small-field inflation models. 

Our general result (2.8) and (2.9) is directly applicable also to non-minimally coupled 
fields, self-interacting fields, and scenarios such as two-field inflation, resonant preheating 
after inflation. One may even generalize the gravity part by not renormalizing to the 
Einstein-Hilbert action, but to different values of a, (3, 62 and include a cosmological 
constant A. 

The work presented here should be seen as an attempt at embedding the dynamics 
of the inflaton and other fields in FRW Universes in a consistent quantum field theory 
framework. It is particularly suited for when the dynamics of the mean field is the main 
concern, but the two point function can also be computed. For many applications (such as 
non-gaussianity in the CMB and warm inflation), it is necessary to go to higher order in a 
1PI expansion to capture the appropriate physics [45]. One may also "resum logarithms" 
into a renormalization group improvement [12]. Alternatively, one could use the 2PI effec- 
tive action [11, 13, 14], but as we argued in the Introduction, the numerical effort at NLO 
and the complications with renormalization currently make this option less appealing. 

There are unresolved questions to do with including also 0(R 2 ) (four derivatives) in the 
effective Friedmann equations [29], allowing for symmetry breaking potentials, Tr? a , < 0, 
and including odd-power terms in the potential acj), cr 3 , (j) 3 , (j)cr 2 , acj) 2 . In connection with 
inflation and (p)reheating, issues related to defining particle numbers relative to some 
specific vacuum persist. It is likely that both the Bunch-Davies vacuum in strict de Sitter 
and the adiabatic vacuum in FRW [20] need to come into play, but a smooth transition 
from one to the other needs to be addressed. 

In this context it may also be important to consider the effect of ignoring boundary 
terms. Although spatial boundaries can safely be sent to infinity, we are dealing with an 
initial value problem, and the initial condition corresponds to applying a localized source. 
This has been seen in the 2PI framework to lead to additional singularities (see [46] for a 
way to deal with this). This calls for closer scrutiny. 

From a technical point of view it would be of interest to relate the renormalization 
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prescription used here to adiabatic regularization of the semi-classical Priedmann equation. 
The truncation at order R 2 is reminiscent of the truncation at four-derivative order in the 
adiabatic regularization, and certainly the fact that all divergencies are cancelled originates 
in this correspondence. The heat kernel approach provides a straightforward way to write 
quantum correction directly in terms invariant operators and their counterterms only. This 
does not follow quite as simply at the level of the equations of motion through adiabatic 
regularization, although presumably we have used one of a number of prescriptions for 
dealing with the finite parts of the counterterms appearing in that context (see also [47] ) . 

Also going beyond R 2 or even including the term 03 in the expansion of the heat kernel 
in a discussion of convergence. Finally, it would be interesting to compare the heat kernel 
result to what one would get from the Keldysh-Schwinger approach in terms of perturbation 
theory. 

The huge success of inflation and the high precision anticipated from current and 
planned CMB measurement mission motivates a more detailed computation of the primor- 
dial perturbations and the dynamics of quantum fields in the early Universe. This requires 
us to go beyond free-field theory and properly include interactions in the dynamics of the 
"classical" fields, and perform correct renormalization. We believe that complementary 
to the powerful but also more cumbersome 2PI formalism, a very useful and intuitive ap- 
proach to this is in terms of the 1PI effective action and variational equations derived from 
it. 
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A Heat kernel method 

This section follows closely [5, 48]. In order to derive a solution for the heat kernel used 
in section 4 we first use and exponential formula to write the variation of the one-loop 
contribution (4.4) as 



dr 5G~ e 



-1„-%tG- 



dr 



Tre 



tG- 



■ (A.l) 



This implies the definition for the heat kernel 



-Id xJ — q I — e 

o / T 



-irG- 1 



dT 



- I dTxy/=g / — K(t;x,x). (A.2) 
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From the exponential form of the heat kernel the following properties can be deduced 

d 

i— K(t; x, x 1 ) = G _1 iT(r; x, x) (A.3) 

OT 

and 

lim K(t; x, x') = 5( ~ x ( A .4) 
r^o ^-g 

Motivated by the flat space definition for the delta function we write the ansatz for heat 
kernel (in n dimensions) 

K(t;x,x') = 1 A l '\x, x'y^'y^nir- x, x'), (A.5) 

where o~(x,x') is the geodetic interval or world function and A 1 / 2 (x,x') is the Van Vleck- 
Morette determinant which can be shown to satisfy 

a^af = 2a, (Acr ; % = nA, [a] = 0, [A 1 / 2 ] = 1, [a ilt ] = 0, 

[A 1 / 2 ;^] = 0, [a.^] = 9lxu , [A 1 ' 2 ,^] = -R„ u , [a^p] = 0, 

r D 2 A i/2i _ nR . R 2 . RgfosK*^* _ R a pR af3 (A g s 

L J 5 36 30 30 ' { ' ' 

where we used the notations 

f(x,x') = f, lim J{x,x')= [/], V„/(s, a/) = /.,,. (A.7) 
The condition (A.4) implies that 

fi(r;a:,a/) = H-O(r), (A.8) 

so if we write O as a power series 

CI(t;x,x') = ^^a k {x,x')(iT) k , (A. 9) 

fc=o 

we kan use (A.3) to write a recursion formula for all the coefficients a&. However, as was 
introduced in [49] and proven in [35] for an operator of the form G^ 1 = — □ + X we can 
sum all the R and X contributions by writing our ansatz (A.5) as 

K(t;x,x') = 1 A^ix, x'y^'^-^'^iT- x, x'), (A.10) 

where 

(M 2 )' = X(*')-^, (A.11) 
o 
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where has also a power series expansion of the form (A. 9). Inserting (A. 10) into (A. 3) 
we get the following recursion relation for the coefficients 



<rfa kitt + ka k = A- 1 /2(A 1 /2 afe _ 1 ) ; ^ + 

With the help of the above formula and the relations in (A. 6) we can solve for the first 
three coefficients 

[oq] = 1, [ai] = 

[a 2 ] =-ha(x-^\ + ^(uR + R a p lS R a ^ s -R a pR^\ (A.13) 

Using this method beyond the first few orders becomes increasingly cumbersome and other 
more efficient methods have been devised [50] and [51]. It is important to realize that 
at no point of our calculation did we choose a particular vacuum in which to define our 
propagator. This means that the result for (A. 10) is not sensitive to boundary condition 
choices, which may in some instances be significant. 



R' 



+ (X' - X) 



(A.12) 
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B Results for counter terms 



For the renormalization conditions in (4.30) we get the following results: 
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where log(/2 ) is defined in (4.14). 
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C Variational formulae 



In this appendix we state the needed formulae for calculating variations with respect to 
the metric. They can be derived with the following identities 



V = -g v<T g w bg pv , 
l 



6^9 



-gg^v Sg 



[IV 



5T1 



g 



^7^ = ^(7=^). 



(C.l) 
(C.2) 
(C.3) 
(C.4) 
(C.5) 



We also make use of the fact that the fields vanish at infinity and that the covariant 
derivative satisfies the Leibniz rule. The variations of the geometric tensors up to mass 
dimension four are 

1 8 f 1 

J d i x^Rf(x)= [--Rg^ + R^-VpVv + g^nlfix), (C.6) 

7=^7 / d 4 x^ R 2 f(x) = [ - l -R 9ilv + 2iV " 2V M V, + 2g^D)Rf(x), (C.7) 



d A x^~g RuvR^fix 



/(*) 



R a 0R a0 g„u + 2R tip Rpj(x) 



+9 lu N p V s (BP 6 f{x)) - 2VPV v (R m f(x)) + n(R pu f(x)), (C.8) 



7=^7 / d^x^g- R a ^ s R aa7 sf(x) = -^R a ^ S R a ^sg^ + 2R/ a ° R upaa f(x) 



+4V rT V p (R fiaup f(x)). 



(C.9) 



If in (C.6), (C.7), (C.8) and (C.9) we choose f(x) = 1, we can use the Bianchi identites 
and commutator formulae for the covariant derivative to get the standard results [6] 



G 



P v 



-gbg^ I d R = ~2 Rg » u + Rpv ' 



(CIO) 



(i) 



H. 



S f 1 

= J d 4 x^g~ R 2 = --R 2 g pu + 2R^R-2V^V u R + 2g pu nR, (C.ll) 
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^H. u = -^==— — — f <fx^ R^Rau 

= --RapR^g^v + IRpuwR" 1 - V,V M i? + ^DRg^ + DR^, (C.12) 

and 

(C.13) 

In four dimensions the higher order tensors are connected via the Gauss-Bonnet theorem 

1 ' y d 4 xV=g (R^R^as + R 2 ~ ^ v R^u) = 0, (C.14) 

fl^ = - {1) H^ + 4^^. (C.15) 



D Geometric tensors in FRW 

When we choose our metric to be 

g^dx^dx" = -dt 2 + a(t) 2 dx 2 , (D.l) 

i.e. of the Friedmann-Robertson- Walker type, we will make frequent use of the following 
tensors 

/2 = 6 f°j + °Y C 2 = 0, G = 24^, (D.2) 
\a z a J a 6 



i 2 

a \ a 2 a 



R Q0 = -3~, ^ = a 2 ( 2— + - ) , (D.3) 



V Vo + sooD)/(t) = 3^ Q /(i), (D.4) 
- V;V; + <feD)/(i) = - fl2 ( 2 ^o + d 2 )/(t), (D.5) 



(1) #oo = -54^ + 36^ -18^ + 36^^, (D.6) 
a 4 cr or a z 



. i , n / d 4 d 2 d a 2 da® a' 4 ' \ 

= a 2 - 18-j + 72^ - 18-j - 24^- - 12 , D.7 

a 4 a 6 a z ar a 



( 2 >iZoo = -18-T + 12^-- 6-5 + 12—=-, D.8 

a 4 a d a" 1 a^ 

f9i t t 9 / « 4 , a 2 d a 2 da® , a® \ _ , 

= a 2 - 6^ + 24^- - 6-5 - 8—=- - 4 . D.9 

\ a 4 a- 1 a z a^ a / 

Having the Gauss-Bonnet theorem (C.15) means that we will not have any need for an 
expression for H^ v . 
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